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Abstract 

The Teukolsky Master Equation is the basic tool for study of perturbations of the Kerr metric in 
linear approximation. It admits separation of variables, thus yielding the Teukolsky Radial Equation 
and the Teukolsky Angular Equation. We present here a unified description of all classes of exact 
solutions to these equations in terms of the confluent Heun functions. Large classes of new exact 
solutions are found and classified with respect to their characteristic properties. Special attention is 
paid to the polynomial solutions which are singular ones and introduce coUimated one-way-running 
waves. It is shown that a proper linear combination of such solutions can present bounded one-way- 
running waves. This type of waves may be suitable as models of the observed astrophysical jets. 

1 Introduction 

At present the study of different types of perturbations of the gravitational field of black holes, neutron 
stars and other compact astrophysical objects is a very active field for analytical, numerical, experimental 
and astrophysical research. Ongoing and nearest future experiments based on perturbative and/or nu- 
merical analysis of relativistic gravitational dynamics are expected to provide critical tests of the existing 
theories of gravity [I] . 

The study of perturbations of rotating relativistic objects in Einstein GR was pioneered by Teukol- 
sky [2] by making use of the famous Teukolsky Master Equation (TME). It describes the perturbations 
sVE'(i, r, 0, lyj) of all physically interesting spin-weights s = 0, ±1/2, ±1, ±3/2, ±2 to the Kerr background 
metric in terms of the corresponding Newman-Penrose scalars. The pairs of spin-weights s with oppo- 
site signs a = sign(s) = ±1 correspond to two different perturbations with opposite helicity and spin 
\s\ — 0, 1/2, 1, 3/2, or 2. Under proper boundary conditions for the TME one obtains quasi-normal modes 
(QNM) of the Kerr black holes. Various significant results and additional references can be found in [3]-[B]. 

The key feature of the TME is that in the Boyer-Lindquist coordinates one can separate the variables 
using the ansatz 'i/{t,r,9,(p) = e^*"*e™'''S'(0)i?(r), i.e. looking for solutions in a specific factorized 
form. Thus, one obtains a pair of two connected ordinary differential equations for the nontrivial factors 
sSuj,E,m{&) and sRLo,E,m{r) - the Teukolsky angular equation (TAE) [2l[3j[7] 

ii^^ (sin0-^,S'„^£,„,(0)^ + .W„,£,m(^)sS'c.,B,m(0)= 0, (1.1a) 
sWu:,E,7n{d) ^ E + Q^iJ- cos^ Q - 2saujcos0~ (m^ + 2mscos9)/ sin^ 9; (1-lb) 
and the Teukolsky radial equation (TRE) [1 13] 



A-^^ (^'^'^ sRu:..E,m{'r?j + .K.,£,™(r) .i?.,£,™(r) = 0, (1.2a) 

, , 1 , 1 dA , , 

sV^.E.m{r) = -rK^-ts- — K-L. 1.2b 
A /A ar 
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The azimuthal number may have values m — 0, ±1,±2, ... for integer spin, oi m — ±1/2, ±3/2, . . . 
for half-integer spin [4]. In Eq. (II. 2[) we use the expressions A = — 2Mr + a^, K = uj{r^ +a^) — ma, 
L — E—s{s+l)+a'^u!'^—2mauj—Aisujr. The real parameter a = J/M is related with the angular momentum 
J of the Kerr metric, M being its Keplerian mass. The two complex parameters uj and E ~ the constants 
of separation, are to be determined using the boundary conditions of the problem. 

The negativity of the imaginary part w/ = 9(a;) < of the complex frequency uj = ujr + iuj ensures 
linear stability of the solutions in the exterior domain of the Kerr metric with respect to the future time 
direction t — >■ +oo [HIS]. In the interior domain the solutions to the TME are not stable [S]. 

From a mathematical point of view the function s3Ct^,£;,m(i, 9, (p) ^ g-'"tg"nv gSu,,E,m{0) sRuj.E.mir) 
actually defines a factorized kernel of the general integral representation for the solutions to the TME: 

s^{t,r,9,^)^^ jdLo jdEY,Au.^E,^e~''^'e''^^,S^^EA^)sRu,E,ra{r). (1-3) 

The formal mathematical representation (|1.3p is written ad hoc as the most general superposition of 
all particular solutions. In it a summation on all admissible values of the two separation constants lu and 
E is assumed. Its usefulness will be illustrated by different examples in what follows. 

It is well known [TU] that the Carter separation constant (which is equivalent to the constant E, 
used here) may be related with the total angular momentum of the solution s'ii{t,r,9,ip). Under proper 
boundary conditions for the TAE this momentum has discrete values defined by an (half) integer I [2]. 
If we are interested in superpositions s'^{t,r,9^ip) of solutions with a definite total angular momentum, 
the integration with respect to the constant E must be replaced with summation over the (half)integer 
Thus, instead of the most general linear mixture (jl.3p we have to use the representation of the solutions 

,*(i,r,0,^) = i- /'d^^^.^„,/,™e-^"*e''"'^,^„,z,™(0),i?^,z,™(r), (1.4) 

introduced in the problem at hand for the first time in ^ . The transition from the representation (|1.3p to 
representation (jl.Sp is formally equivalent to the use of a singular kernel proportional to the sum of Dirac 
(5-functions: 6(E — sE{uj, I, m)) in (II. 3p . Here sE{uj, /, m) belongs to some spectrum which is specific 
for the given problem and is completely defined by the corresponding boundary conditions, see section 7. 

Further on, the boundary conditions may fix some discrete spectrum for the frequencies uj in (jl.Sp . 
Then the integral on uj will be replaced by discrete summation over some a;„. This is equivalent to the 
use once more of a singular kernel, now proportional to S{uj — ujn). As a result one obtains 

,*(t, r,9,^)=Y,Y.Y. e""" sSn,l,n^{9) sRnJ,m{^). (1-5) 

n I m 

The lack of a rigorous mathematical theory explaining the Teukolsky separation is an old and well- 
known problenfl. The study of the QNM not only illustrates the above situation but also shows 
that the kernel s^u,E,m{t,r,9^ip) can be singular with respect to the variable r at infinity and at the 
horizons. In the existing literature only regular with respect to the variable 9 kernels s^uj,E,m{t,r,9,ip) 
are in use. In the present paper, we start the consideration of both regular and singular with respect 
to the angle 9 kernels in the integral representation (|1.3I) of the solutions to the TME. Different types 
of kernels are to be used for solution of different boundary problems. Note that from a physical point 
of view the regularity of the very solution s'^{t,r,9,Lp) in equations (|1.3I) . (11.41) and (|1.5p is important. 
The kernels like s^ij,E.mit, r, 9, ip) are auxiliary mathematical objects. One is often forced to use singular 
kernels in the natural integral representations of the solution to physical problems. The regularity of the 
very physical solution s^(i, f, 9, ip) with respect to the variable 9 depends on the choice of the amplitudes 
sJ^uj,E,m- It can be guaranteed by a suitable choice of these amplitudes, as shown in section 10. 

Despite the essential progress both in the numerical study 11 of the solutions to equations (|l.lap and 
(ll.2ap and in the investigation of their analytical properties |il2 , at present there exists a number of basic 

^The author is grateful to unknown referee for this important remark. 
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questions remaining unanswered. For example, it has been well known for a long time |13) that the TAE 
(|l.la|) and TRE (|1.2al) can be reduced to the confluent Heun ordinary differential equation p3] written 
here in the following simplest uniform shape |171 118] : 

H"+(a + ^ + 2±l]H'+(^ + ^]H = 0. (1.6) 

\ z z— 1/ \z z— 1/ 

The constants /i and v in Eq. (jl.6p are related with the constants a, /3, 7, (5, 77, accepted in the notation 
HeunC(a, /3, 7, S, 77, z) as follows: 

/3 + 7 + 2 

= /i + — a 2 1 (1.7a) 

= 2 2 ■ ^ ' 

To the best of our knowledge we still do not have a detailed description of the exact analytical solutions 
to the TAE (II. lap and TRE (|1.2ap in terms of the confluent Heun function HeunC(Q!, /?, 7, 5, 77, z) - the 
unique particular local solution of Eq. (|1.6p which is regular in the vicinity of the regular singular point 
z = and obeys the normalization condition HeunC(a, ^, 7, 5, r;, 0) = 1 [T4]|l. Note that other particular 
solutions to equation (|1.6p , as well as its general solution, are not termed " confluent Heun's functions" , 
according to the accepted modern terminology [13]. The reason is that, in general, other solutions can be 
represented in a nontrivial way in terms of solutions HeunC(a, /?, 7, S, f], z) of the corresponding arguments. 
Hence, from a computational point of view it is sufficient to study only the Taylor series of this standard 
local solution and its analytical continuation in the complex plane Cz. Thus, the instrumental use of the 
confluent Heun function HeunC(a, /?, 7, ^, 77, z) - the basic purpose of the present paper, is much more 
advantageous than the simple fact, recognized already in [13], that the TRE and TAE can be reduced to 
the confluent Heun equation (11.61) . 

In the late 2006 a program for filling the above gaps in the study of the TME was started as a natural 
extension of the papers [15], where a similar approach was developed for the Regge- Wheeler equation 
(RWE). The first results were quite stimulating [16], but serious difficulties came across in both analytical 
and numerical studies. This is because the theory of Heun's functions, as well as numerical tools for 
calculations with them still are not developed enough. 

Here we pay special attention to the polynomial solutions of Eq. (|1.6p . According to [M], the confluent 
Heun function HeunC(a, /?, 7, (5, 7;, z) reduces to a polynomial of degree iV > of the variable z, if and 
only if the following two conditions are satisfied: 

- + ^^+N+l = 0, (1.8a) 
a 2 

Ajv+i(^) = 0. (1.8b) 

We call the first condition (|1.8aP a "(5Ar-condition" , and the second one (ll.8bp - a " AAr+i-condition" . An 
explicit form of the " A^v+i-condition" in form of a determinant useful for practical calculations, as well 
as a novel derivation of confiuent Heun's polynomials can be found in [17j . A recurrent procedure for 
calculation of A^r+i (|1.8bp and its relation with Starobinsky's constants are presented in |18j . 

On the other hand, the so-called algebraically special solutions to the RWE and the TRE were dis- 
covered long time ago [19]. These are of a generalized polynomial type, i.e. products of polynomials and 
simple non-polynomial factors which are elementary functions. According to the existing literature, these 
solutions describe pure incoming or pure outgoing waves. The algebraically special solutions still are not 
discussed in terms of Heun's polynomials. To our knowledge, attempts for application of this class of 
solutions to real physical problems cannot be found in the existing literature on gravitational physics. 



^ In the present paper, we use the Maple-computer-package notation for the Heun functions. Basically, this notation is 
borrowed from the two mile-stone papers on modern theory of Heun's functions by Decarreau et al. in 1141 : at present it 
seems to be most popular, since the Maple package is the only one for analytical and numerical work with Heun's functions. 
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The only exception are the recent papers [161 120] . where one can find some preliminary results. There, 
special polynomial solutions of the TAE were considered in more detail. In particular, in the first two of 
the articles [50] it was demonstrated that these singular with respect to the angular variable 9 solutions 
can describe collimated waves which resemble the observed astrophysical jets. In the third of these papers 
the spectrum of electromagnetic jets from Kerr black holes and naked singularities in the Teukolsky 
perturbation theory was calculated for the first time, using some of the basic results of the present paper. 

Very recently the algebraically special solutions of the RWE and the TME were proved to be relevant for 
the study of instabilities of different kind of some more or less "exotic" solutions to the Einstein equations 
|21) . Physical manifestation of the instabilities of the mathematical solutions are the explosions of the 
corresponding objects. Therefore, it seems natural to look for a perturbative description of explosions in 
terms of solutions of the TME, which are stable in the future and instable in the past. The confluent 
Heun functions give a rigorous mathematical basis for analysis of these problems. 

On the other hand, the recently found properties of confluent Heun's function [17] show that one 
can introduce a new subclass of " JAr-confluent-Heun's- functions" , which obey only the (5Ar-condition - 
Eq. (I1.8a|) . In [18 is shown that such "(5Ar-solutions" of the TRE and the TAE define the most general 
class of solutions, for which properly generalized Teukolsky-Starobinsky's identities exist. Moreover, this 
approach reveals the existence of Teukolsky-Starobinsky's type of identities for Regge- Wheeler and Zirilli 
equations, as well. Here we study in more detail the (5jv-solutions to the TRE and TAE. In particular, 
we show that the regular solutions to the TAE, which are the only class of solutions to the TAE, used up 
to now [UI31[3) are precisely nonpolynomial ^Ac-solutions. In contrast, the polynomial solutions to the 
TAE of all spins are shown to be singular around one of the poles {9 ~ 0, or 9 ~ n) oi the unit sphere 

(2) 

Sg ^ and regular around the other one. This new situation reflects the specific properties of the confluent 
Heun function. It is not consistent with our experience, based on the work with hypergeometric functions, 
solving the angular part of the Laplace equation in celestial and quantum mechanics, or in electrostatics. 
It is well known that in the last case solutions regular in the interval 9 E [0, it] (including both the poles) 
are polynomial. 

In the limit a — > 0, when the Kerr metric approaches the non-rotating Schwarzschild one, there exist a 
smooth transition from perturbations of the Kerr metric to perturbations of the Schwarzschild metric in 
terms of the Weyl scalars, but a simple transition from the solutions of the TME to the solutions of the 
RWE (see [IS1[TS]) is not possible [3]. Nevertheless, the mathematical analogy between the corresponding 
solutions becomes quite transparent when the solutions are represented in terms of the confluent Heun 
functions [TS] [T5]. The limit a — is traced in more detail in section 4.1.2 - for the TRE and in section 
7 - for the TAE. 

This way, using confluent Heun's function, we hope to obtain a more clear picture of the quite compli- 
cated present-day state of the arts in the perturbation theory under consideration and its possible further 
developments and new physical applications. 

The main purposes of the present paper is to report some of the basic results, obtained for a detailed de- 
scription of the exact solutions of the TME in terms of the confluent Heun function HeunC(Q;, /3, 7, (5, 77, z), 
to introduce a large number of new classes of such solutions, and to formulate some interesting boundary 
problems for the TRE, TAE, and TME in terms of confluent Heun's functions. 

Besides the already stressed new developments, in the present paper for the flrst time we introduce 
and study the differential invariants of the Weyl tensor, which indicate in an invariant way both the 
event and Cauchy horizons of the Kerr metric as singular points of the TRE (section 2.1), the explicit 
form of 16 classes of exact solutions to the TRE in terms of confluent Heun's functions (section 2.2), a 
new classiflcation of the solution to the TRE, based on specific properties of confluent Heun's functions 
(section 3), especially, the class of (5jv-radial solutions and, in particular, two unknown infinite classes of 
exact solutions with equidistant complex spectra of frequencies, two novel classes of polynomial solutions 
to the TRE (section 4), the explicit form of 16 classes of exact solutions to the TAE in terms of confluent 
Heun's functions (section 5), a new concomitant confluent Heun's function and its application to the 
TAE (section 5), a new classification of the solution to the TAE, based on specific properties of confluent 
Heun's functions (section 6), especially, the class of (^Ar-angular solutions, a novel description of the regular 
solutions to the TAE in terms of confluent Heun's functions (section 7), two classes of singular polynomial 
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solutions to the TAE (section 8), 256 classes of exact factorized solutions to the TME (section 9), an 
explicit construction of exact bounded solutions to the TME with spin 1/2, using the singular kernel, 
built from the polynomial solutions to the TAE (section 10), and novel general exact solutions of the 
TME in the form of one-way running waves (section 10). It seems natural from physical point of view to 
use these solutions for study of the relativistic jets. Some other general conclusions and perspectives for 
further developments are outlined in the concluding section 11. 

2 Exact Solutions to the Teukolsky Radial Equation in Terms 
of the Confluent Heun functions 

2.1 Explicit form of the TRE and Geometrical Character of its Singularities 

Much like in the case of the Schwarzschild solution, for the Kerr one we have a complicated space-time 
structure and a different physical meaning of the space-time coordinates in the different domains. For 
example, consider, as usual, only the real values of r. In the interior of the Kerr metric: < r_ < r < r_|_ 
- between the zeros r± = M zL V — a? < M of the function A (i.e., between the Cauchy horizon 
r_ and the event horizon r+), two of the eigenvalues: At and Aj. of the metric in the Boycr-Lindquist 
coordinates simultaneously change their signs. Indeed, one pair of eigenvalues is \e = gee = r'^ -\- a? cos^ 9 
and \r = grr = {r"^ + a^cos^9)/A. The second pair of eigenvalues is the roots At, A^ of the equation 
•^^ ^ (fftt + 9<t><t>)^ + 9tt94>4> — gt4,- Their product equals AtA^ = — Asin^ 6. The last expression, together 
with the form of grr proves the simultaneous change of the signs of the two eigenvalues At, A^, when the 
variable r crosses the horizons r±, since the determinant g = — (r^ -I- cos^ 9)^ sin^ 9 of the metric does 
not vanish there. As a result, between the two horizons r± the variable Un — x (z (— oo, oo) plays the role 
of the interior time and the variable ri„ — t is the interior radial variable. We use the following Kerr- 
metric-tortoise-coordinate: = r -I- a_|_ In |(r — r+)/(r+ — r_)| — a_ In |(r — r_)/(r+ — r_)| e (— cx),cx)), 
where a± = r±. It is a straightforward generalization of the tortoise variable for the exterior domain 

r S (?"+, cxi) proposed in [2]. Since our expression is valid in the interior domains, too, the inverse function 
defines r = r{ti„) when r S (r_,r_|_). In the second interior domain r <t- the variables r and t restore 
their original meaning. For a detailed analysis of the light cones in the Kerr geometry see [22]. This 
consideration is necessary for understanding of the physical meaning of the solutions to the TME in the 
different Kerr-space-time domains. 



The explicit form of the TRE 



+ (1 + s) 



dr^ \ r — r+ r — j dr 

1 

/I 1 \ , 

w ia^ + — am 



\ 



(r — r-i-) (r — r_) \r — r — r. 



-E + s{s + l)-a^uj''' + 2mauj + 'iisujr\ ^^^^ =0 (2.1) 

(r — 7'_|_)(r — ?■_) 

shows that it has three singular points: r ~ r± and r — oo. In the present paper, we consider only the 
non-extremal Kerr metric with real r^ > r^ > 0. Then the first two are regular singular points, and the 
third one (the physical infinity r — oo) is an irregular singular point. The symmetry of Eq. p.l[) under 
the interchange r+ ^ r_ is obvious. Thus, we see that the two horizons of the Kerr metric are singularities 
for the TRE which arc to be treated on equal footing. Do these singularities have an invariant meaning 
independent of the coordinate choice? 

It is well known that the algebraic invariants of the Riemann curvature tensor 3lijki are not able to 
indicate the horizons of the Kerr black hole and one usually considers them as pure coordinate singularities 
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of the metric in the Boyer-Lindquist coordinates. In contrast, the circle r = 0,9 = it /2 is a, singularity of the 
algebraic invariants of the Riemann tensor [3]. Since the pure algebraic invariants of the tensor Jiijki do not 
fix completely the geometry, their consideration is not sufficient to recover all gauge-invariant space-time 
properties. For this purpose one must consider a large enough number of high-order-differential-invariants 
of the Riemann tensor ^5] . 

It is not difficult to find differential invariants of the Riemann tensor of the Kerr metric which are able 
to distinguish both the horizons r± and the ergo-surface gu = 0. Indeed, let us consider the following 
algebraic invariants of the Weyl tensor Wijki'. h = ^"Wyfc; W*-''""' - the density of the Euler characteristic 
class, and I2 = j^'Wijki - the density of the Chern-Pontryagin characteristic class [M]. Let us put 

/ \l/6 / \l/6 

ih - ihY'"^ = A = |A|exp(iV)- Then r = f^j cos(V'/6) and p ^ i^A cos(?/;/6)'i are obviously 
invariants of the Weyl tensor - nonalgebraic and nondifferential ones. In the Boyer-Lindquist coordinates 
one obtains p = r -I- — cos 6 and gu = I — 2M/ p. The differential invariants of first order 



D/i = -(Vlnr)^ = — I 1 - — + ^ ) , (2.2a) 
rp \ r 



1 ( 2M 



Dh = (V In pf - (V In rf = 1 - l) (^1 - (2.2b) 

indicate the two Kerr black hole horizons, the ergo-surface and some other geometrical objects in the 
Kerr space-time. Thus, the two horizons r± of the Kerr metric are shown to be invariant objects, being 
singularities of the same kind in equation (|2.ip . 

In the limit a — >■ we have p — >■ r and the differential invariant in equation (j2.2bp becomes trivial: 
DI2 — > 0. In the same limit, the differential invariant (|2.2ap produces a nontrivial result (1 — 2M/r) jr^ 
for the Schwarzschild metric which is similar to the one derived already in the third of the papers |23| , as 
well as in the very recent seventh one. Differential invariants similar to (2.2) are considered in the sixth 
of the references [23] without any application. 

2.2 Explicit Form of the Local Solutions to the TRE 

The analytical study of the solutions to the TRE and TAE was started in [5^ and continued by different 
approximate methods [5J[T2 without utilizing of Heun's functions. Using the confluent Heun function one 
can write down 16 exact local Frobenius type solutions to the TRE (|2.1I) in the form: 



«<£,m,.„,.„a,('^;'^+,^-)A^/'=e^-^0^^^z^-^HeunC(a„a^,a^/3^,a^7±,5±,^/±,^J, (2-3) 

which is very similar to the form of the solutions to the RWE [T^l US] • Her^ 

a+ ^ sau,E,ra{r+,r^)^2iu{r+ - r_)= i2puj/Q.a, (2.4a) 

/3+ =sPu,,E,m{r+,r^)= s + i{m-uj/n^) Ip, (2.4b) 

7+ =slu,E,m{r+,rJ)= s-i{7n-uj/Q.+) /p, (2-4c) 

<5+ =s5u,E,ra[r+,r^) = a+ (s - iuj{r+ + r_))= a+ {s - iuj/Qg) , (2-4d) 

2 2 2TO2f]2_^2 (2mf7a-a;)2 if ajfl+V , , 

T — r_ T-\- — T — A 

z+ = z+{r;r+,r^) = , z_ =:z_(r; r+, r_) = ; + = z+z_ = - (2.4f) 

r+ — r_ r+— r_ (r+— r_)^ 



^Note that the notation z± in Eq. I l2.4fl l is consistent with the Hmits z± — > ±00 for r — > 00. Their relation with the 
notation of the parameters of the Kerr metric r± is illustrated by the equations z±{r^\rj^,r^) = 0. The labels it in the 
notation in Eq. I I2.3I I are related with the labels of their arguments z±, not with the labels of the parameters rj-. 



6 



The discrete parameters aa,<Jp,a-y have values ±10. In equations (I2.4al) - (l2.4el) we use the following 
quantities: the angular velocity of the event horizon fl-^ = a/2Afr+ = \Jr^ jr^ j (r+ + r„), the angular 
velocity of the Cauchy horizon — a/2Mr^ — VM-/*^/ (^+ + the arithmetically-averaged angular 
velocity Vta — (ri+ + 57_)/2 = l/(2a), the geometrically-averaged angular velocity Vtg = ^JVL^SlZ ~ 
1/(2M), and the new dimensionless parameter 

p=i (V^W^- V^~a7) - ^/n^)^^AP/a^-i. (2.5) 

Note that the inverse relation r± = y^l^~/^l±/ -I- permits us to replace r± with r2± wherever it 
is necessary, thus making transparent the duality of the parameters r± and r2±, as well as the behavior of 
the above quantities under interchange of the two horizons: r-|- ^ Q,^ ^ p i— > — p, i^a.g ^a,g 

- invariant. 

The parameters a_ , /3_ , 7_ , (5_ , •q_ can be obtained by interchanging the places of the two horizons: 
r+ ^ r_ in p.4ap - p.4ep . This procedure may be substantiated using the known properties of the 
confluent Heun function under changes of parameters [14j . One can check directly that this way we 
obtain indeed solutions of equation (j2.ip . 

According to equations (|2.3p and ()2.4fj) . the behavior of the solutions e m a a (^5 around 

the corresponding singular points z — ^± (t*^ ; r+ , ) = is defined by the dominant factor (z^)'^'^'^*^^. 
All other factors in equation ()2.3p are regular around these points. The same solutions are in general 
singular around the corresponding singular points z = z± (r± ; r_|_ , r_ ) = 1. 

Only two of the sixteen solutions (|2.3p are linearly independent. Nevertheless, it is necessary to know 
all of them since for different purposes one has to use different pairs of independent local solutions. 

Using the known asymptotic expansion of the confluent Heun function [14 we obtain two asymptotic 
solutions of Tome type. These are local solutions of the TRE around its irregular singular point \r\ = oo 
in the complex plane C^: 

- e--K(^-+-)'"0 ^a, , ao = 1. (2.6) 

The notation ±cx) in (j2.6p denotes the two directions: r — !■ +oo and r — > — oo on the real r-axis for 
approaching the irregular singular point \r\ = oo in the complex plane C^. For the coefficients aj = 
o-j,u,E,m,aa,tTf,,a^ One has a recurrence relation |14) which shows that they increase together with the 
integer j. Hence, the asymptotic series (|2.6p is a divergent one. 

As seen from (EH), sR~^E,m,a^,a^,a-,ir;r+,r^) = si?+^£;,,„,^„,^^^<,^ (r; r_ , r+). Hence, one can introduce 
a new parity property of the solutions and construct symmetric and anti-symmetric (with respect to the 
interchange r+ r_) solutions of the TRE: 

sRu^E^.in.a„.a^,a-,{r;r+,r^) = i (^sRXE,m,a^,a^,a^i'r'''r+^r^) + sRZ.E,m^a^.a^,a^ir-,r+:r-)) , (2.7) 
sR^^E,rn,aa,<yf,.a^{''"T'''+T'''-) ~ ^ (''-^i-E^m.fo ,(7^ ^+ ' ) ^ s-^u;,£;,m,(T„ ,(T^ ,cr^. (^i ' )) ■ 

Clearly, these solutions are singular at both horizons in the general case. When one considers the two- 
singular-point boundary problem |14l on the interval [r_, r+] in the Kerr black hole interior, the solutions 
(|2.7p may be regular at one, or at both the ends for some values of the separation constants cj and E. Since 
this boundary problem is still not studied, at present we are not able to make more definite statements 
about this case. 

^Further on, tTj, =sign(x) denotes the sign of the real quantity x. The only exception is a = Us where we skip the index s. 
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3 A New Classification of the Solutions to the TRE, Based on 
the (^TV-Condition. Novel Radial ^^v-Solutions 



For the TRE the (5Ar-condition reads: 
where 

We call radial 6 n -solutions the solutions defined via the (5jv-condition p.ip . 

The calculation of the values of the coefficients in equation p.ip yields two very different cases: 

1. In the first case ±±— -^±±± = ^ind we see that one is not able to fix the frequencies ± ± ± 
and ± ± ±- Instead, choosing CTq = cr^ = cr^ = —a and using (|3.ip one fixes the non- negative integer 
N in the form 

,A^ + 1 = 2|s| > 1 for |s| > 1/2. (3.2) 
Thus, the degree of the polynomial Ajv+i-condition is fixed, too. 

2. In the second case the coefficients are nonzero and one can fix the values of the frequencies 
sLd^ g.^ g.^ from equation (|3.ip . Thus, one obtains two different types of exact equidistant spectra: 

a) For ±,± = ^^,±,± = ±2, M+ = ^t,±-± = ^ and ?^|,±,± = ^^,±,± = (iV + 1) the 
(^Ar-condition (13. ip produces the pure imaginary equidistant frequencies with iV > ~ integer: 

s^N,m,T,±.± = ^''N,m,T,±.± = ±* + ^) / = ±i ^(^(iV + l)/2. (3.3) 

Note that these frequencies depend neither on the spin-weight s and azimuthal number m, nor on the 
rotation parameter a. The spectrum is not influenced by the rotation of the waves and of the very Kerr 
metric. The frequencies p.3p are defined only by the monopole term in multipole expansion of the metric. 

b) For all other cases the coefficients „ „ , and „ „ are not fixed and one obtains the fol- 
lowing two similar double-equidistant spectra of frequencies with > - integer and m - (half)integer: 



s^7V,m,T,T,± " ^^iV.m.T.i.T ^ ^+ ± ip{N + 1 =F s)) ; (3.4a) 
^^W,m,±,T,± = «^W,m,±,±,T = ^- ± ip{N + 1 ± s)) . (3.4b) 

A set of important new mathematical properties of the radial (5Ar-solutions can be found in jl71 118] . 
In [T^ it is shown that these solutions define the most general class of solutions to the TRE for which 
the properly generalized Teukolsky-Starobinsky identities exist. The solutions which satisfy the relation 
(|3.2p were studied in [2l [3] without utilizing the Heun functions and the ^Ar-condition. The last condition 
turns to be valid automatically for the solutions to the TRE studied in [2l [3] . The infinite series of the 
solutions with equidistant spectra p.3p and p.4p are introduced and considered for the first time in the 
present paper. In the third of the papers |20j one can find an interesting and unexpected recent application 
of the formulas p.4p for fitting of the spectra of electromagnetic jets from Kerr black holes and necked 
singularities. 



4 Polynomial Solutions to the TRE 

The (^AT-condition yields the basic classification of the solutions described in the previous section 3. As 
a result, one obtains two classes of polynomial solutions to the TRE, imposing in addition the A^r-)-!- 
condition (jl.8bp . In what follows we will use the determinant form of the A^v+i-condition given in 17J. 
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4.1 The First Class of Polynomial Solutions to the TRE: 



The solutions of this class correspond to the first case in section 3 and obey equation (|3.2p . The inequality 
sN = 2|s| — 1 > excludes the existence of scalar perturbations (|s| = 0) of the first polynomial class. 



4.1.1 The General Case: 

For brevity, we denote the solutions £; „ _c, -o- -cr(^' _e m(''' ^-)- ^O'" them the 
parameter /x takes the values /i = ^/i^ j, ?'-), fc = 1, . . . , 2|s| - the solutions of the algebraic equation 

(ll.Sbp . which now takes the form: (^) = 0. Its degree is 2|s| = 1, 2, 3, or 4, depending on the spin of 
the perturbations \s\ = 1/2, 1,3/2, 2. Making use of (|1.7bp . and (j2.4ap - (|2.4ep . we obtain for the separation 
constant £' = si?^ J, ^^j^(r+, r_), A: = l,...,2|s| the expressions 

sE^M,Jr+,r^)^,fi^, ,^ir+,r^) + \s\i\s\ - 1) - au;{au; - 2m) + 2ia(2|s| - l)^r^, (4.1) 

Applying the explicit expressions for the roots sM^j, m(''+'^-)' ^® obtain: 

,S±^(r+,r_)=-aV + 2au;m-i : for |s| = i, m = ±1/2, ±3/2, . . . ; (4.2) 



(r+,r_)= -aV + 2aw (m - (-1)'' ^1 - m/au?j : for fc = 1, 2, |s| = 1, m = ±1, ±2, . . . (4.3) 



For the gravitational waves (|s| = 2) one has to find the quantities sA*^ ^ ?'-) solving algebraic 

equation of the fourth degree A^(/x) = 0. The explicit form of its roots is too complicated and not 
necessary for the purposes of the present paper. It is more instructive to demonstrate here the result, 
obtained using the Taylor series expansion of the solutions s/^^ ^ j„(^+, ^-) around the zero frequency 
a; = 0. 

Thus, we obtain for \s\ = 2, fc = 1, 2, and m — ±2, ±3, . . . the eight series of values: 



/ 3A'I\ 



m■^+^i-lr2m ( (^^-1)^ + ^) (aa.)^+ (4.4) 



+03(aw). 



For \s\ —2, k — 3,4, and m — ±2, ±3, ... we have other eight series of values: 



^ = i {-l)''4y/mauj(l + i3 ^1 + (^^^ - mauj +©2(0^) j + 



+8maw - 6 ( 1 + ( ^ g j ^ j +03(aw). (4.5) 

Clearly, these series describe two kinds of solutions with a completely different behavior around the 
origin a; = 0. In particular, the series (|4.4p and (|4.5p have different limits: 2 and 0, respectively, when 
a; — )■ 0. For the solutions (14. 5p the origin w = is a branching point, etc. 

The independence of the values of sE'^j^^n (|4.4p and (|4.5I) on the upper labels (±) is a result 
of the polynomial character of the solutions, i.e. of the regularity of the corresponding HeunC-factor 
simultaneously on both the horizons r±. 

For a complete solution of the problem one has to determine the frequency u. Hence, one needs an 
additional relation between the parameters E and uj. This relation may appear when one solves the TAE 
(See the next sections 5-8.). 

The first class of polynomial solutions to the TRE is introduced and studied in detail for the first time 
in the present paper. 
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4.1.2 The Special Case of the Schwarzschild metric: 



For the special value of the parameter a — we have r_ — 0, — 2M. This is the case of perturbations 
to the nonrotating Schwarzschild black hole described in terms of the Weyl scalars. For simplicity, here 
we use units in which 2M — 1. The parameters in the solution (I2.3P acquire the limiting values 

a+ — 2iuj, /3+ = s, 74. = s + 2iw, (5+ = 2iuj{s — iuj), 77+ = —E + — ; 

a_ — — 2ia;, /3_ = s + 2iLu, 7_ ~ s,S^ = —2iLu(s — iu!),i]^ = —E + + 2w + 2isLjj. (4-6) 

These differ from the values of the parameters of confluent Heun's functions in the Regge- Wheeler approach 
to the perturbations of the Schwarzschild metric |15| . 

In the limit a — >■ equation p.ip does not define the frequency w, if (Ja = T^/? = ±cr-y = — cr, because 
then one obtains L^^j^^^^ = 0. If, in addition, a — sign(s), then the JAr-condition is fulfilled for the 
special polynomial solutions of the first class denoted as ^i?^ ^ ^^^(r) = ^ _^ ^cri^\ I7 0)- Equation 
p.ip yields the relation gN = \s\ — 1 > 0. Scalar perturbations of this type do not exist. 

In the case of integer spins \s\ — 1,2 the roots /i = s^J■^ k m' ^ = 1, • • • , |s| of the equations (/i) ~ 0, 
(|1.7bp . and (|2.4ap - (|2.4ep with r+ = 1, r_ = and a — produce the following simple expressions for 

= : for \s\ = 1, (4.7) 
sE^ f^ „^ = 1 - (-1) Vl - i6(7t^ : for \s\ = 2, fc = 1, 2. (4.8) 

For a complete solution of the problem, one needs an additional relation between the parameters E and 
ui. This relation may be found by solving the TAE, see sections 5-8. 

The above considerations of the limit a — >■ and the corresponding results for the Schwarzschild black 
hole in terms of confluent Heun's functions are new and obtained for the first time in the present paper. 



4.2 Second Class of Polynomial Solutions to the TRE: 

According to the results of section 3, the solutions of this class originate from the second case of the Sn- 
condition and fall into two subclasses: a) and b). The complete definite frequencies ^oj^ ^ ^ ^ - formu- 
lae (13. 3p and p.4p . yield algebraic equations A^^^(^) = with (A^+1) roots fi = „ m cr„ o-^ a-, ('''+''''-)' 
n = 0, 1, . . . , A^. It seems difficult to derive explicit analytic expressions for these roots, but their numer- 
ical values can be easily obtained. Using the values of sfJ'^nmac af, o-^ (''-1- > ) ^^"^ equations (jl.7bp . 
(|2.4ap - (|2.4ep we obtain complete definite values for the parameter E = sE^ nmrjc rrn a., (^+: ''-)• 

a) In the case of frequencies p.3p we obtain 

sE^,^,m..,^a.-. = s\J^N,r..ra,a,^.,-. + I'^KI-'^I " 1) + «^(3a^ " 2m) + 4 + 2*ac. (2Af | s | - r±). (4.9) 

b) In the case of frequencies (I3.4ap . (j3.4bp we have, respectively: 

uj^ r\.\ I -ip^ ^ (4.10a) 

^<„.™,+,+,--.A*^,„,m,+.+,- = .^^^.„,™.-,-,+ -.A*^,„,™,-,-,+ = ±*2(m+2aa.(l-2AfVa2))/2p- 

-i2psauj - A[m^ -f 2m (l-3MVa^) auj - {l-bU'^/a^) {aujf^ / Ap^ . (4.10b) 

With Lu and E given by equations (13.31) . (13. 4p and (14.91) . (I4.10p we have no more free parameters in 
the problem at hand. As a result, the corresponding solutions to the TAE are fixed unambiguously by 
the designated group of equations obtained for the second class of polynomial solutions to the TRE. This 
situation is completely new, unexpected and described here for the first time. 
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5 Exact Solutions to the Teukolsky Angular Equation in Terms 
of the Confluent Heun functions 

In terms of the variable x = cos9 the TAE has three singular points. Two of them: x- = —1 (i.e., 9s ~ tt 
- South (S-)pole) and x+ = 1 (i.e., 6'jv = - North (N-)pole) are regular singular points. The third one 
Xoc = c» is an irregular singular point. It is remarkable that introducing the notation 

z+ = z+(6l) = (cos(6l/2))^ z- ^ z-{e) ^ {sm{e/2)f , z+ + z_ = 1; (5.1) 



and 

a± = zLAauj, b± = s ^ m, c± = s ± m, d± = zLAsauj, n± — =p 2sauj — a^ui'^ — E. (5.2) 

we can write down 16 local solutions of the TAE in the form 

a_j_ z c_|_ 

sS'^,E,m,a^,a^.a^ = e'"" Z^^" — — HcuuC (cTaOi , Cr&6± , CTcC^ , d± , , Z± ) (5.3) 



which is very similar to the form of Eq. (12. 3p . 

Following the corresponding properties of the TAE (II. lap [2] , the solutions (|5.3p have the symmetries 

-s'S'^,_E,m,cr„,cri,,cr^(^ ^ ^) ^ ,E ,m~ <y ^ -at, - a ^ (5.4a) 

Note that according to Eq. (|5.3p . the behavior of the solutions sS[^ ^ around the correspond- 

ing singular points z = z+{9s) = Z-{6n) = is defined by the dominant factor (2;^)°^''''*^^. All other 
factors in (15.31) are regular around these points. The same solutions are in general singular around the 
corresponding singular points z = z^{9n) = z^{9s) = 1- Hence, at this point we have a complete analogy 
with the case of the TRE. 

Only two of the sixteen solutions (j5.3l) are linearly independent. Nevertheless, it is important to know 
all of them, since for various purposes one can use different pairs of independent local solutions, see below. 
If one chooses some two linearly independent solutions, then one can represent the other fourteen using 
this basis. Unfortunately, at present the form of the corresponding coefficients is completely unknown. 

We can establish simple relations between some of the different solutions ()5.3p in proper domains of 
the parameters s and m, if we divide the whole plane {s, m} into four sectors. In each of them we choose 
the solutions with the same regular asymptotic behavior around the corresponding pole as follows: 

I. Sector s > 0, \m\ < \s\: 

s^ui,E,m(^) ~ s'5'^,_E,m, +++ = s'5'tj,_E,m, — \-+~ '^^(^,E,m,++- — sS^^^^j^^_^_ g"^^ ['^'^^ 2 I ,(5.5a) 



II. Sector m < 0, \s\ < \m\: 

( 9\^^^™ 

s^u],E,mi^) — sS^Ejm, ~ s'^i.j,_E,m,H —sS^^.m, h ~ « "^cj^B^rn.H ^ 6)2^0 I ''^^ 2 j •(5-6b) 

III. Sector s < 0, \m\ < \s\: 

s^ij,E,mi^) — sS^^E„^^ =s'5'„_£;_„i^-| = i'^w.B.m, '>^'^,E,m,^ ^ e^jr \ 2 j 7(5-7a) 



^^u],E%ii^) — sS^^^ ,^^ —sSuj,E,rn,-{ " s'^w,_E,m, h^'"^w,S,m,H h g'^Q \^^^ 2 ' •(5.7b) 
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IV. Sector m > 0, \s\ < \m\: 

s^uj,E,mi^) — sS^^^ ,^^ =sS^^^ j^ _^ ^sS^^^ j^^^ h^'''^tJ,-B,™,H ^ g"^^ 2 J j(5-8a) 

/ 0\ s+m 

s^uj,E,mi^) — 'S'5'w,£;,™.,+++ ^s'S'^^B^m^ — \-+ ^ '^^'^ ,E ,m, + + - ^ sS^^^ „^ __^__ g2^Q ( ^^'^ 2 j ' 

Note that in each sector the four solutions in the above relations of type (a), or in the above relations 
of type (b) are equal, since under standard normalization the local regular solution around any regular 
singular point of the TAE is unique. 

In the case of the TAE there exist an additional complication. The numbers s and m are simultaneously 
integers, or half-integers. Then /3 — abb_^ — CTb(s=pm) is an integer and, in particular, it may be a negative 
integer. However, the confluent Heun functions IIeunC(a, /3, 7, (5, ry, z) are not defined when /3 is a negative 
integer |14j . Therefore, if /3 = erf, 6^ < is a negative integer, we must write down the corresponding 
solutions in the form 

sS^,E,,n,a^,a„a, ^ ^I' ^ ^ HcunC (cr^ fl^ , cTfefe^ , a^c^ ,d^,n^,z^). (5.9) 
For this purpose we define the concomitant confluent Heun functiorj^ 

HeunC fg. /S, 7, S, t], z) = z-^HeunC(a, -fi, 7, 5, ry, z) / ^ ^-dC- (5.10) 

J (HeunC(a,-/3,7,<5,r?,C)) 

This function is well defined for negative integer /3 — a^h^ < 0, together with the confiuent function 
IIcunC(a, — /?, 7, 5, ry, z). In this case, the function z^'^HeunC(Q;, — /?, 7, 5, 77, z) represents the local regular 
solution to the confluent Heun equation (11.61) around the singular point z = and the concomitant 
confluent function HeunC(a, /3, 7, (5, r/, z) represents a second linearly independent local solution, which is 
singular around this point. We need the concomitant Heun function to construct second independent local 
solution in the case /? = 0, too, since in this case z~^HeunC(a, — /3, 7, J, 77, z) = HeunC(a, /3, 7, 5, 77, z) can 
not be used for this purpose. 

It can be shown that for negative integer /3 the concomitant confluent Heun function has the form 

l/3| 

IIeunC(a,^,7,5,?7,z) = V — + /7i(z) + /i2(2)ln(z), aU c„ 7^ 0. (5.11) 

z — i 2" 

n=l 

Here hi 2{z) denote two definite functions of the complex variable z which are analytic in the vicinity of 
the point z — Q. In the problem at hand = |/3±| = |s =F m\. The logarithmic term is present in the 
concomitant confiuent Heun function when — 0, too, but then we have no poles in the solution (I5.11[) . 
For = its form otherwise is similar to ()5.1ip . One can reach the last results using general analytical 
methods described, for example, in [25]. 

The above consideration sows that in the case of the TAE we can construct only eight local solutions 
(|5.3p which are single-valued functions of the variable z. The other eight solutions, being in the form 
(|5.9p . are infinitely- valued, because of the logarithmic term in Eq. (I5.1ip . 

6 A New Classification of the Solutions to the TAE Based on 
the ^TV-Condition. Novel ^^y- Angular Solutions 

For solutions r? „ „ „ (15.31) to the TAE the (5jv-condition reads: 

= =Fm hiV+l+ (CTaH Is. (6.1) 

^Notc that for any value of the parameter /3, when the confiuent Heun function in the right hand side of Eq. 115.1011 is 
well defined, its left hand side represents a second, linearly independent solution of the confluent Heun equation. 
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We call angular SN-solutions the solutions defined via the (5Ar-condition (|6.ip . To some extent these 
solutions are similar to the radial (5Ar-solutions introduced in section 3. A set of important new mathemat- 
ical properties of the angular i^Ar-solutions can be found in [T71 [T^ . In [TH] it is shown that these solutions 
define the most general class of solutions to the TAE for which properly generalized Teukolsky-Starobinsky 
identities exist. 

Comparing equation (|6.ip with the corresponding one for the TRE - p.ip . we see both essential 
differences and similarities. For the coefficients in equation (|6.ip . which are analogous to the ones in (|3.ip . 
one obtains: 



S. 



Hence: 

i) The coefficients 0-5 vanish identically, in contrast to the coefficients in equation (|3.ip . 
Consequently, there are no cases in which the condition (j6.ip can fix the frequencies uj. 

ii) The form of the coefficients of both equations (|3.ip and (|6.ip is the same only for M ja = -\/2- 

iii) The coefficients H'^ of both equations are of the same form. 

We obtain two different cases depending on the coefficient (db — CTc) in front of the azimuthal number 

m: 

1. The first class angular (5Ar-solutions with CTc = o^, and Oa = <yh = <^c = ^f- As a result. Eg. ()6.ip 
fixes the degree of the second polynomial condition A^r+i = in the same form as equation p.2pl : 

siV + 1 = 2|s| > 1 for |s| > 1/2. (6.2) 

In the case of the TAE the set of (5jv-solutions consists of ones with integer parameters /? — <7\,hj_ 
of both signs. According to Eqs. (j5.9p and (|5.11l) . for negative integer cr;,6^ < in the solutions we 
have logarithmic terms. Such solutions are infinitely- valued functions. To exclude this physically not 
admissible case, one must impose the additional requirement (7(,6^ > 0. As a result one obtains <J<Jm = ±1, 
(Tftb^ = |m| — |s| > and OcC^ = —\m\ — \s\ < 0. Thus, the single-valued angular Jat- solutions of the first 
class with spin |s| > 1/2 correspond to sectors II (|5.6p and IV (|5.8p and acquire the form 



s<i.,™(^J=e-'^"'^"^±(^±) ' (^,) ' X (6.3) 
xIIcunC(— lo-maw, \m\ — \s\, —\m\ — |s|, 4(Tmaw|s|, 2(Jmau)\s\ — a^uP — E, z^). 



Obviously, the solutions (j6.3p are regular around the singular points =0. Their behavior around the 
singular points z^ = 1(<^4> z^ — 0) is more complicated. We can study this behavior using the expansion of 
the solutions ()6.3p with respect to the basis of the two linearly independent local solutions ()5.9p around 
the points =0, which are well defined for aam = ±1, o'bb^ — \m\ — |s| > and <7cC^ — — \m\ — \s\ < 0: 



S^,E,mi^J = ,r±(^,i?,m)e2-— . (z^)^- (z^)^ X (6.4) 



X HeunC (4am auJ , —\m\ — \s\, \m\ — \s\, —4(7mau}\s\, h 2(T„iacj|s| — a uj — E, z^) + 

|.r.| + |s| 

xHeunC(4crmaw, |r7i| -I- \s\, \m\ — \s\, —4:amauj\s\, h 2(T„jaa;|s| — a^cj^ — E, z^). 

Now it is clear that in the general case, when sr^(u;, E, m) ^0, the solutions (|6.3p are singular around 
the corresponding points z^ — I and in addition - still infinite valued, because of the poles and of the 
logarithmic terms in the concomitant confluent Heun function in Eq. (|6.4p . as well as because of the 



"The alternative case (T[, = Cc = —fa leads to a uninteresting relation A'' + 1 = 0. 
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singular factor (2^) ^ ■ One can remove at once all these unwanted properties from solutions 
imposing the condition (w, E', to) = 0. Unfortunately, the explicit form of the connection constants 
sTf2i^j to) is completely unknown. At present, this is one of the main unsolved problems in the theory 
of the confluent Heun functions. 

Another way to avoid the logarithmic terms in the solutions (|6.3p . ()6.4p is to impose the A^r+i- 
condition, reducing this way confluent Heun's functions to polynomials. We consider in detail these two 
possibilities in the next sections 7 and 8. 

2. The second class angular J^v-solutions: ai, = —Oc- Then we obtain 

sA^m,o-„,o-6,-crt, + 1 = ± mCTb - OaS > 1. (6.5) 

Now the additional requirement (7b{s T ™) ^ and ()6.5|) yield the solutions 

.<s,™,-.,.,-.(^J =6^^— ± (zj^ (zj"^ X (6.6) 
xHeunC(=F4(Taa;, |s| =F ctto, — |s| =F ctto, ±4(Taa;|s|, ^ =p 2craLLi|s| — a^w^ — i?, z^). 

with |s| > 1/2 and to restricted in the asymmetric finite intervals 1 — |s| < icrm < |s|, i.e., — |s| + (1 ± 
(T(Tm)/2 < |to| < |s| — (1 =F o'(Tm)/2. These solutions have ctq = — ct^ (Tc = — cr, sNm,~a-,cr.~a + 1 = 
|s| ± crm > 1 and correspond to pairs {s,m} in sectors I (|5.5p and IV (15. 7p . Under the above conditions 
the solutions sS^ ^ „^_^ ^ _^{z^) (|6.6I) are obviously regular around the points = 0. Their behavior 
around the second regular singular points — 1 can be studied using the following expansion with 
respect to the corresponding local basis 

sS^^E..m,-..a.-J^J = sTt{u;,E,m,-a,a,-a)e^'''^-^^ (zj"^ (zj^ x (6.7) 



xHeunC(±4(Taw, —\s\ =F am, \s\ =F am, ^4:aauj\s\, ± 2craa;|s| — a^u)^ — E,z^) 



^rf{oj,E,m,-a,a,-a)e^'^'''^^^ (z^)' {zj 



i2, ,2 
X iu 



o2 



xHeunC(±4(TaCi;, \s\ ± am, \s\ =F am, =p4(Taw|s|, ± 2(Taa;|s| — a uj — E, z^). 

Since — 2|s| < —\s\ =F am < — 1 and 1 < |s| ± am < 2\s\, the two independent local solutions in (|6.7p are 
well defined. The first solution (with concomitant Heun's function) is singular around z^ = 1 (<^4> z^ = 0) 
and, in addition, infinite-valued. One can remove at once all these unwanted properties from solutions 
(|6.7p imposing the condition sr^(a;, _E, m, — ct, ct, — o") =0. Unfortunately, at present the explicit form of 
the connection constants sr^2('^i E, to, —a, a, —a) is completely unknown, too. 

Another way to avoid the logarithmic terms in the solutions (|6.6|) . (j6.7p is to impose the A^r+i- 
condition, reducing this way confluent Heun's functions to polynomials. We consider in detail these two 
possibilities in the next sections 7 and 8. 

As seen, in the case of the TAE the only role of the (5Ar-condition is to relate the degree N of the 
Ajv_|-i-condition with the spin-weight s and the azimuthal number to and to select the proper solutions. 

Note that up to now only regular solutions to the TAE, which obey the condition (j6.2p have been 
studied and used in the literature [2l[3l[7]. In section 7 we develop a new approach to the regular solutions, 
based on confluent Heun's functions. The nonregular angular ^Ar-solutions, subject to the condition (|6.2p . 
and the infinite series of solutions, subject to the condition (|6.5p . are introduced and considered for the 
first time in the present paper. 

7 Regular solutions of the TAE 

The spectral conditions s^f{^,E,m) = and s^f{Lo,E,m,—a,a,^a) — ensure the regularity of the 
solutions (|6.3p and (|6.6I) , as seen from the formulas (|6.4I) and (|6.7I) . One is not able to use these conditions 
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directly, since the explicit form of the connection constants (a;,i<^,m) and (w, i?, m, — cr, ct, — cr) is 
not known. Therefore, we are forced to use a roundabout way to find the regular solutions to the TAE. 
Suppose we have a solution sS^^^^{9) which is regular around the S-pole {6s = tt) and another 

-^{0) which is regular around the N-pole (6'jv = 0). We will have a solution s'S'J'f '^„(6'), 



solution ,5* 

regular everywhere in the interval G [0,7r], if and only if m(^) " const x sS~^^^^{9), i.e., if the 



Wronskian vanishes: W 



0. This condition determines the constant E in the 
form E = E{auj, s, m, I), I being a (half)integer. The Wronskian will vanish for any 6 G [0, tt], if it is zero 
for some Oq € (0, tt). 

To utilize this idea for all values of the parameters s and m, we have to divide the whole plane {s, m} 
into four sectors and to choose the solutions sS^^^^j^{9) defined by Eqs. (|5.5p - (j5.8p . 

The spectral condition makes equal the solutions of group (a) and the solutions of group (b) in each 
sector. It can be written in different equivalent forms combining in pairs one solution from the group (a) 
and another one from the group (b). Below we give the simplest form of this condition in each sector, 
written here for the first time in terms of confiuent Heun's function HeunC and its derivative HeunC'. 
The set of all conditions (|7.ip defines the separation constant in the whole plane {s, m} in the form 
E = [vn? + s^)/2 — a^w^ + e{auj, m, s). The new parameter £(acj, m, s) is to be found from the following 
transcendental equations: 
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(7.1a) 



= - in sector I, 



(7.1b) 



- in sector II, 



(7.1c) 



- in sector HI, 



(7.1d) 



= - in sector IV, 



valid simultaneously for all values of £ (0,7r). Thus, the two-singular-points boundary problem for the 
TAE is solved. It yields a countable set of values E{auj,m,s,l) numbered by some (half)integer I: I is 
integer for an integer spin, or half-integer - for a half-integer spin. Due to the symmetries (j5.4l) of the 
solutions to the TAE, the different relations (j7.1ap and (I7.1cl) . or (j7.1bp and (j7.1dp give similar results. 
More precisely E(auj, to, — s, /) = E{aLu, to, s, I) and E{—aLu, — m, s, /) = E(auj, m, s, I). 
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An important consequence is that all the regular solutions obtained this way are angular (5Ar-solutions 
with the same (|6.2p in sectors II and IV, or with the same sA^m aa.a-b (|6.5[) - in sectors I and III. This is 



because between the solutions (|5.6p and (|5.8|) we certainly have ^Ar-solutions: sS^ ^ „j , for s > 0, and 

s'^^,£;,m.+++' fo'' * < 0- Between the solutions (|5.5p and (|5.7p (5Ar-solutions are sS'J^; -p_|_-p, for m > 0, 
and sS^ ^ ^ T±T' ^'^'^ JTi < 0. As a result of uniqueness of the regular solutions with given values of the 
parameters, all regular solutions inherit the i5Ar-property. Hence, all regular solutions of the TAE obey 
the Teukolsky-Starobinsly identities [18]. 

Let us consider the limit aw — > of the regular solutions to the TAE. Since 

HeunC(0, /3, 7, 0, r/, z) = (i_2)/3+7+i+V^^+7^+i-4», ^ 

/^ /?+7+l+V/3^+7^ + l-47? /3+7+l-V/^^+7^ + l-4ry \ 

I 2 ' 2 ;p+J^;^ I: I'-^j 

in this limit the Heun functions in Eqs. (|5.5p - (|5.8p and (|7.ip can be reduced to the Gauss hypergeometric 
ones. Then, using the well-known properties of the Gauss hypergeometric function 2F1 one can derive from 
Eqs. (|7.ip with aw = the spectrum E{0, s, Z, m) = + I = l{s, m, I) = max(|m|, |s|) + ?^, [ = 0, 1, 2, . . . 
The values of the separation constant E{0, s,l,m) in this case are real. The numerical analysis of Eqs. 
(j7.ip written directly in terms of confluent Heun's functions confirms this standard result for the limit 
auj = 0. The corresponding regular confluent Heun functions in sS^^qi ^{9) in this case are reduced to 
Jacobi's polynomials - in the case of an integer spin, and to their spin-weighted generalizations - for a 
half-integer spin 

The solutions E{au}, s,l,m) for small auj and integer spin have been studied many times in the 

form of Taylor's series expansion E{auj, s,l,m) = ^(^ + 1) + X]^i Ej^s,i,m{0'Ujy without use of Eqs. (|7.ip 
and without utilizing the Heun functions. A little bit surprising thing is that the solutions sS^i'^{9) 
with auj ^ 0, regular at both poles, are not polynomial and can be represented as an infinite series with 
respect to Jacobi's polynomials. Here we describe the regular solutions to the TAE in terms of confluent 
(5Ar-Heun's functions for the first time. 



8 Polynomial Solutions of the TAE 

8.1 Singularities of the Polynomial solutions to the TAE 

The polynomial solutions to the TAE are a special subclass of the angular (^Ar-solutions studied in section 
6, since both of the two conditions (|1.8p are valid for them. Being a polynomial in z, the HeunC-factor is 
regular at both regular singular points = 0, tt. Then the singularities of the polynomial solutions around 

the poles are defined completely by the factors (z^)'^'"''*^^ and (z^)'^"'^*^" in Eq. ()5.3p . Thus: 

1. In case of the first class angular (5Ar-solutions (|6.3p with 

|m| > \s\ 

we see that the singularities are defined by the factor (z^)"^^"^*^^ which gives sS'^ e m^'^±) ^ 

(sin|) i_e., singularity at the N-pole 9 = 0, and sS~ ^ mi'^±) ^ ('^o^f) i.e., singu- 

larity at the S-pole 9 = tt. 

2. In case of the second class angular (5Ar-solutions (16. 6p with 

-|s| + (1 ± aa^)/2 < \m\ < \s\ - (It (Ta^)/2 

we see that the singularities are defined by the factor (z^)'^"^*^^ which gives sS^ ^ .^^ _^ ^ _^{z^ ) ^ 

(sin|) (1'*!+'^™^^ jg^^ singularity at the N-pole 9 = 0, and sS~ ^ _^ ^ _^{z^) ^ (cos |) i.e., 
singularity at the S-pole 9 = tt. 
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As a result, we see that in any case the polynomial solutions are regular around one of the poles and 
singular around the other one. 

Using relations (|1.7b[) and (|5.2I) we obtain the general formula for the constant E in the form 

~fi^ — auj^ ^ 2(Ta(lT + {<7a + (Tb)s)auj H — -m (ufcrn T 1) H — -s (crbS+ 1) • (8-1) 

Further analysis shows that some of the properties of the two classes of polynomial solutions to the 
TAE resemble the corresponding properties of the two classes of polynomial solutions of the TRA, but 
there exist also some essential differences. 



8.2 First Class of Polynomial Solutions to the TAE: 

These are the solutions sS^ ^ ^ _^ _^ with AAr+i-condition fulfilled. For them the specific requirement 
abb^ > yields the restriction |m| > \s\ and the condition ()6.2|) is fulfilled independently of the values 
of the azimuthal number m. As in the case of the first class polynomial solutions to the TRA - section 
4, the value s = is eliminated by (|6.2p . Hence, we have an infinite series of the first class polynomial 
solutions to the TAE for all admissible values of s and m. Preserving the style accepted in the previous 
sections we denote the polynomial solutions to the TAE of the first class as sS[^ ^ ^ = sS[^ ^ „ _^ _^ 

For them the A7v+i-condition reads A2[s|(/x) = and has 2 |s|-in- number solutions slJ-^km- From 
formulae (|8.1|) one obtains 

= s^^t,k,m + - 1) - a^{auj - 2m) T 2a(2|s| - l)auj, (8.2) 

where fc = 1, . . . ,2|s|, s = ±1/2, ±1, ±3/2, ±2 and in addition |m| > \s\. 
Solving the Ajv+i-condition, we obtain for the different values of \s\: 

sE^„^^ -a'^w'^ + 2aujm-^: for |s| = i, m = ±1/2, ±3/2, . . . ; (8.3) 



+ 2auj (^m - (-l)''v/l - m/au?j : for fc = 1, 2; |s| = 1, to = ±1, ±2, . . . (8.4) 



The values (|8.3|) and (|8.4p of the separation constant E obtained for the first class polynomial solutions 
to the TAE are the same as the corresponding values ()4.2p and (|4.3p for the first class polynomial solutions 
to the TRE. Important consequences of this unexpected fact are considered in a separate paper |26) . 

For the gravitational waves (|s| = 2) the quantities fc m solutions of the algebraic equations of 
the fourth degree AJ(/z) = 0. We do not need here the exact form of these roots. It is quite complicated. 
Below we present only the form of the separation constant E for the TAE obtained making use of the 
Taylor series expansions of the roots around the point auj — 0. 

Thus, we obtain for \s\ = 2, fc = 1, 2, and to = ±2, ±3, . . . the following eight series of values: 



»E^ ] 



- Amauj ~ i{~lfl2^{m - l)m{m + 1) {aujf^ + 6 [rn^ - {aujf + 0z/2{auj), (8.5) 
and for |s| = 2, fc = 3, 4, and m = ±2, ±3, . . . another eight series of values: 
s£^^_fe_„ = -(-l)''4VTOaw ^1+ ^3to- ^ au + 02{auj)^ +?>mauj-Q{m^ - ^{aujf -VOziauj). (8.6) 

As seen, for gravitational waves of the first polynomial class the values (j8.5|) and (j8.6|) of the corre- 
sponding constants E differ substantially from the analogous values (|4.4p and (|4.5I) of the constants E 
obtained for the TRE in section 4.1.1. This is in sharp contrast to the case of neutrino waves (|s| = 1/2) 
of the first polynomial class and to the case of electromagnetic waves (|s| — 1) of this kind. 

It can be shown that this phenomenon refiects the difference between the Starobinsky constants for 
solutions with spin 2 to the TAE and for solutions with the same spin 2 to the TRE p] |3l [18]. The 
solutions to the TAE and to the TRE with the same spin 1/2 or 1 have the same Starobinsky constants. 

Despite the above essential difference, the first-polynomial-class-solutions to the TAE and to the TRE 
with spin 2 have similar qualitative properties, discussed at the end of section 4.1.1. 
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8.3 Second Class of Polynomial Solutions to the TAE: 

We have a finite number of second class polynomial solutions to the TAE for which the relation ac ^ —<Jb 
holds. For brevity, we list here only the ones of spin 2, 1 and 1/2. For them the conditions > and 
<7bb± >0 must be satisfied simultaneously, yielding the requirement — |s| < — |s| + (1 ± aam)/'2 < \m\ < 
\s\ — (1 =F crcrm)/2 < |s| - almost opposite to the analogous requirement |m| > \s\ for the polynomial 
solutions of the first class. Altogether there exist only the following 32 polynomial solutions of the second 
class sS^ ^ m =F ± T ^^^^ ®P^^ ^' ^ ^'^'^ ^1'^' 

s = +2, m = -l, 0, 1, 2; s = +l, TO= 0, 1; s= +1/2, m= +1/2, -1, 2; 
s = -2, m = -2,-l, 0, 1; s = -l, m = -l,0; s= -1/2, m= -1/2, +1/2; 
s = +2, m = -2,-l, 0, 1; s = +l, m = -l, 0; s= +1/2, to= +1/2,-1/2; 
s = -2, m = -l, 0, 1, 2; s==-l, TO== 0, 1, s= -1/2, m= -1/2, +1/2. (8.7) 

The relation between the constants E and lo follows from (18. ip . when fi in it is replaced by the solutions 
of the Ajv_|_i-condition in the form A^_|_^|(/i) = 0. Here we omit these relations. 

9 The 256 Classes of Exact Factorized Solutions 
to the Teukolsky Master Equation 

Combining solutions to the TRE and to the TAE studied in the previous sections we can construct the 
following 256 classes of exact factorized solutions to the TME 

s^tt^.,.^..^,,.^-,,.^...^...^^^^^^^^^^ (9-1) 

For specific physical problems one has to impose specific additional conditions, like stability conditions, 
boundary conditions, casuality conditions, specific fixing of the in-out properties, regularity conditions etc. 
Thus, one selects some specific combinations of solutions to the TRE and to the TAE in Eq. (|9.1|) and 
derives the spectrum of the separation constants u) and E in the given problem. 

For example, choosing solutions to the TRE which enter both the event horizon and the 3D-space 
infinity we study the Kerr black holes (for a < M), or necked singularities (for a > M) [11 [3]. If in 
addition we choose regular solution to the TAE, we will obtain the standard QNM of the Kerr black 
holes, or necked singularities. Using in Eq. (|9.ip other solutions to the TRE and/or to the TAE, we may 
hopefully describe different physical objects and phenomena, for example, collimatcd jets, see in [Tallin]- 

The solutions (|9.1I) do not necessarily have a direct physical meaning, see the Introduction. Instead, 
some linear combination of the specific solutions, which obey proper boundary conditions, is to describe 
the Nature. In general the solutions (|9.1|) have to be considered as auxiliary mathematical objects - 
(maybe singular) kernels of integral representations (|1.3|) of the physical solutions. The choice of the 
corresponding amplitudes sAuj,E,m,<Ta,<Tf3,a-^,<Ta,a-b,a-a will fix completely the physical model and may ensure 
the convergence of the integrals and discrete sums to physically acceptable solutions. Since at present we 
have no rigorous mathematical treatment of this complicated issue, we will study it in the next section 
using some constructive examples. 

10 Construction of Bounded Linear Combinations of Polynomial 
Solutions to the TAE 

We have seen in section 8 that the polynomial solutions to the TAE are singular and unbounded with 
respect to the angle 9 around the N-pole, or around S-pole. These solutions produce a singular kernel 
in the integral representation (II. 3p . It is important to know whether it is possible to have bounded with 
respect to the angle 6 e [0,7r] solutions s^'(i, r, 6*, (ys) defined by Eq. (|1.3p . despite the singular character 
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of the kernel in it. The answer to this question is a quite nontrivial issue. Here we reach a positive answer 
for perturbations of spin 1/2 in several steps. 

Let us consider the simplest case of double polynomial solutions of the first class to the TME with spin 
1/2 and s = a/2. For them we have an essential simplification, since according to Eqs. (|3.2p and (|6.2p 
sN = 0. Hence, the HeunC-factors in both the radial and the angular polynomial solutions are equal to 
const = 1. The value of the separation constant E = —a^uj'^ + 2aujm — | is uniquely defined in both cases 
by Eqs. (|4.2p and (18. 3p . Hence, the integration over the constant E i n (|1.3p produces only one term with 
this fixed value. As a result, the corresponding singular kernel (|9.ip if|j: 

^•^u, E m(i, r, 9,ip)=6(E + a^u? - 2mauj + 1/4) A'^e""^^" , (10.1) 

^ ' ' VsinS 

and — 0, if ct = +1, or 



where = t + cr (r* - zacos^), Wa = e**'' cot cj)^ = ip + In 



-TT - 



if cr = -1. 



(2) 

The complex variable Wa- defines a stereographic projection of the two-sphere §^ g on the com- 
pactified complex plane Cw„ ■ Its use is critical for further analysis of the problem. Note that after the 
transition from real variables {0, 4>ct} to the complex one one must introduce an additional phase 
factor exp(— zs^cr) in the spin-weighted spheroidal harmonics, due to the back rotation of the basis (See 
the paper by Goldberg et al. in [4].). In the case of spin 1/2 the introduction of such factor exp{'^i(j>cr/2) 
is equivalent to a transition in what follows from half-integer to integer values of the azimuthal number 
m and a replacement m ^ m ± 1/2 in the factor 6 {E + a^u? - 2mau + 1/4) in Eq. (fTOTi) . 

Taking the trivial integral on the variable E, one obtains from the representation (|1.3p and Eq. (|10.ip 

1 



^^{t,r,0,^) = A{r)-^^i\W,\ + \W,\-^)/2 }] | — /dc. e^"^' . A^,™ | (T4^.)" . (10.2) 



m— — oo 



Since in this case we have no other restriction on the frequencies uj, different from the stability requirement 
5(w) < 0, the otherwise arbitrary integration contour G C,^ in (110. 2p must lie in the lower complex 
half-plane. Suppose that the amplitudes ^Ai^^m and the contour L^^ are chosen in such way that for all 
771 G Z there exist well defined integrals 

Jdoj e-''^^' ^^„,„, = ^2l™(r,). (10.3) 

Then 



^■^{t,r,0,^) - A(r)-^V(|W^.| + \Wa\-^)/2 ^ .2l„ (T.) (W^.)" . (10.4) 

?n — — oo 

Suppose, in addition, that in some ring domain jW^j € {\W\' , IM^I"): < \W\' < \W\" < oo the sum 
J2m=-oo f (-^<t) (W^ct)™ = .|2t (To-, Wa) represents a convergent Laurent series of some analytic function 

(Tct, Wa)- For this purpose the coefficients s.^m {Ta) in Eq. (|10.3I) for 777 > and, independently, for 
777. < must satisfy some of the well-known criteria for convergence of the corresponding series. Thus, 
we finally obtain a solution to the TME with spin 1/2 which depends on an arbitrary analytic function 
!L^{Ta,Wa) of the two variables Ta and Wa'- 



^^{t,r,0,if) = Air)- — ^i\Wa\ + \Wa\-')/2 ^%{Ta,Wa). (10.5) 



^To simplify formula l|10.1|l . we have omitted some constant factors in the corresponding solutions to the TRE and TAE, 
which do not depend continuously on the real variables r and 9, but may have different values outside the event horizon, in 
the domain between the event horizon and the Cauchy horizon and inside the Cauchy horizon. This is a legal operation, 
since one can include these factors in the amplitudes s-A... in the representation Ijl.Sp . considering separately the different 
domains where these factors are constant. 
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Returning to the Boyer-Lindquist variables one can check directly that (|10.5I) indeed gives a general 
solution to the TME with spin 1/2. The explicit form of the variable shows that outside the event 
horizon these solutions describe one-way-running waves: outgoing to space infinity running waves - for 
a — —1 and incoming from space infinity running waves - for cr = +1. 

Now it is easy to remove the singularities from the z-axis, i.e. on the poles 9 — O.n. For example, 

let us choose 



has no singularities on the poles 9cj ■ 



2p 



-Wa^) 12. Then !^^{t,r,e,Lp) = /^{r)- ^ / 1 - (f>a sin^ 9^ 
0,7r, but this way we have worked out two new singular lines 



In r-r^ = ±7r/2 on the equatorial plane 9 — Tr/2. Hence, this way the singular line of the 
solution has been only deformed and translated to a new position. The same happens if we choose the 
more general function ^^{T„,Wa) = l/^J {a{T„)W^ + b{T„)Wa^ + c{T„)) /2. In this case, the singular 

z-axis will be deformed, translated and doubled to the non-static singular lines (j)^ = ip - 



2p 



In 



4'i,2 = arg(M^i^2) on the (in general) moving cones 6* = 0i^2 = arctan (|W^i^2|~ ) , where Wi^2 are the two 
roots of the equation a{Ta)Wa-\-b{Tcr)W~^-\-c{Tc) — 0. Here we have chosen a special form of the function 
|.2t (Tcr, Wa) which yields finite nonzero values of the solution on the poles 9 = 0,Tr. 

It is possible to chose the function £. 21 (T^ , Wa- ) with the denominator which is a sum of polynomials 
of higher degree with respect to variables Wa and W~^. Then the solution (|10.5p equals zero at the N 
and S-poles and we can work out an arbitrary number of singular lines of the solution related to the 
zeros of the denominator. At first glance, this possibility may not seem to be interesting for the physical 
applications, since on the singular lines the linear perturbation theory in use is not applicable. We mention 
it here just to have a clear mathematical picture. It is interesting to study the same situation in the whole 
nonlinear theory and to know whether in it the singular lines may be replaced by regular ones. If so, the 
perturbation theory under consideration indicates a possible complicated structure of the exact radiation 
field on the Kerr background. 

The most important question for a correct application of the linear perturbation theory under consid- 
eration, is whether one can find a regular analytical function |.2t (T^, W^) without singularities in the com- 
plex plane Cvk^/{0, c»}, i.e., with the points = and W^, = oo punctured and which, in addition, can 
remove the unbounded increase of the solutions due to the singularities of the factor ^(|Wcr| -f iVKcrh"^) /2 
in (jlO.Sp . We give a positive answer to this question constructing two explicit examples: 



1. Using the basic equality ^ M^™/„j(z) = exp 



for the modified Bessel functions 



Im{z) [57] we choose the coefficients in (|10.4I) in the specific form £.%rn (To- 
where lo — ujn + iuj is a fixed frequency and a — sign(|wij| — |w/|). Then 



^^Ut,r,9,^) = A{r)- — ^i\Wa\ + \Wa\-^) /2 exp(--c^X' 



expli 



(10.6) 



is a stable solution, since by construction it goes to zero when t — > +oo. It is not difficult to obtain its 
limit when — >■ 0, tt in the form 



^Jm^ {^^Ut,r,9,^)) = A(r)-— exp (--c^^W 



X lim 




(10.7) 



Here 



'^ui,a:0,7r = ± Us + — In 

\ 2p 



— a arctan 



t + or* 



+ arg(a;), for 9 = 0, or tt 



(lO.i 



is the hmit of the total phase of the term ^ (Wa-i-W^ ^) wTcr and Ta;o,7T = t + o'(r* =F ia). In Eq. (110.81) the 
sign {+) corresponds to the limit 9cr ^ and the sign (— ) - to the limit 9cr tt. Formula (jl0.7p shows 
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that when T„^o.;o,7r G (^f 7 f ) the solution ^^'^^(t, r, 9, (p) is bounded everywhere in the interval 9 S [0, tt], 
since in this case lim (^.^(^(i, r, 0, (p)) = 0. Otherwise this limit diverges and the solution is singular 

and unbounded around the poles. 

Actually, the value of the parameter Ttj,cr;o,7r is not defined from a geometrical point of view, because 
the value of the angle (p is completely arbitrary on the poles 9cr = 0,7r. As a result, we can choose any 

value of the parameter T^^^-o.tt without changing the geometrical points associated with the N and S-poles 

(2) 

of the sphere §g ^. Since the different values of this parameter yield different solutions of the TAE, we see 
that under the boundary conditions at hand the corresponding differential operator is not self-adjoint |28| . 
but its self-adjoint extensions do exist and can be fixed by suitable fixing of the free parameter T(^^a-;0,7r- 
An analogous phenomenon is well known for the potentials V{x) ^ or in quantum mechanics 

[28]. Note that around the poles 9a — 0,tt the potential in the TAE (jl.ip has precisely the same behavior: 
sWu:,E,m{d) ~ 1/^^ for 9^0, and sWi^^E,7n{d) ^ - 7'')^ for 9 ^ n. In our case, the fixing of the 
parameter Tu:,a;0,7T G (^f i f ) makes the solutions (|10.7p to the TME for spin 1/2 smooth and bounded 
everywhere in the interval 9 G [0,7r], i.e., physically acceptable. 

2. Another solution, which is finite everywhere in the interval 9 Cz [0, tt] but has an infinite number of 
bounded oscillations around the poles = 0, tt can be obtained using the following equality for the Bessel 

fimctions J™(z): ^ (_i)™VK2m (j^(z)f = Jo( {W+W'^) z\ jlT]. Now we choose the coefficients in 

m— — OO ^ ^ 

(Unai in the specific form ^Sla™ (T,) = (-l)'"exp (-fw^Tj^) (J2m(^r,))' and ^^2,n+i (T,) = using 
the same notation as in the previous example. Then 

^vl/„(i,r,0,^) = A(r)-^yaW^TrW^F)72 cxp (-|c^X') Jo{{Wa + W-') luT,) (10.9) 

is a stable solution to the TME with spin 1/2. Taking into account the asymptotic expansion of the Bessel 
function Jo{z) ^ \/ ^ cos(z — 7r/4) we obtain in the limits 9^ — >■ 0, tt: 



\im^^ {^^Jt, r, 9, ip)) = A(r)-^ cxp [-^lu'T,\, 



y/TTioTa:0,Tr 



X lim cos -Uy^ ^ e«T.„;o,, ng.lO 

e„-)-o,7r \ \ sm6' J J 

As seen from Eq. (110. lOp . there exist only two choices of the free parameter: T^^a;0,TT = 0, tt, for which the 
solutions (|10.9I) are finite everywhere in the interval 9 G [0, tt] - a critical property for the use of the linear 
perturbation theory. Approaching these poles the solutions oscillate infinitely many times with bounded 
finite amplitudes. In this sense, the N and S-poles remain singularities of the bounded solutions ()10.9p . 
Moreover, the gradients of the bounded solutions (|10.9p are unbounded around the poles. 

Obviously, superpositions of solutions (|10.6p . or (jl0.9p with different complex parameters w, running 
in some (discrete or continuous) sets in C^j, describe more general bounded solutions to the TAE with 
spin 1/2. 

One more remark. In the case a = +1 the solutions (jlO.Sp are unbounded on the horizons r± due 
to the factor A(r)~^/^. These stationary singularities cannot be removed by any choice of the function 

(Tcr, W^ct), since it depends on the two variables and Wa, not on the single one r. The variables 
t, r, 9 enter in and the variables (p, r, 9 enter in in a complex way. As a result, the variable r cannot 
be disentangled from the the function |.2l (To-, Wo-) and one is not able to compensate the singularity due 
to the factor A(r)~^/^ which does not depend neither on the time t, nor on the angles p and 9. 

11 Conclusion 

In the present paper, we have demonstrated that the confluent Heun functions are an adequate and 
natural tool for a unified description of linear perturbations of gravitational field of the Kerr metric 
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outside the event horizon, as well as in interior domains. These functions give us an effective tool for 
exact mathematical treatment of different boundary problems and the corresponding physical phenomena. 
They can help us to solve old mathematical issues, related to the Teukolsky separation of the variables, 
as well as to study new physical problems. The same approach works, too, for the Regge- Wheeler and 
Zerilli equations in the Schwarzschild metric |15[ IT5] . 

Large classes of exact solutions to the perturbation equations of the Kerr metric were described here 
for the first time. All possible types of solutions were classified uniformly in terms of confluent Heun's 
functions and confluent Heun's polynomials, using their specific properties. As we saw, the variety of the 
different solutions and possible spectra is much reacher than, for example, the variety of the corresponding 
solutions and spectra of the Hydrogen problem in quantum mechanics, solved exactly in terms of the 
confluent hypergeometric functions . Mathematically, this is obviously caused by the presence of one 
more regular singular point in the confluent Heun equation. 

We have to stress especially the newly obtained singular polynomial solutions to the Teukolsky an- 
gular equation. They differ drastically from the well-known analogous regular polynomial solutions to 
the quantum Hydrogen problem. The singular polynomial solutions to the Teukolsky angular equation 
present an auxiliary mathematical construction which seems suitable for simple and natural perturbative 
description of the coUimation of radiated fields of all spins |s| > in the Kerr metric, see [TFl [20] . 

For spin 1/2 we have proved that the singular kernels, constructed from polynomial solutions, can 
produce bounded solutions of the continuous spectrum to the TME with very interesting physical prop- 
erties: These solutions describe collimated one-way running waves in the Teukolsky perturbation theory 
correctly. 

For spin 1 we have also proved the existence of double polynomial solutions of the continuous spectrum 
to the TME and can reach similar results in a more complicated way, since there we meet a new physical 
phenomenon - the electromagnetic superradiance [29l [21 [3] . For spin 2 the supperradiance is known to be 
quite stronger, but we have no continuous spectrum of the TME and the problem needs special treatment. 
We shall consider these two important cases separately. 

One can hope that the collimated one-way running waves, cropping up for the first time in the present 
paper, are able to describe the real astrophysical jets, observed at very different scales in the Universe. This 
still speculative idea needs a more detailed mathematical development and a careful confrontation with 
the real astrophysical observations. It indicates the existence of a new universal mechanism for collimation 
of radiation of all spins |s| > by the pure gravitational field of rotating compact astrophysical objects 
of different nature. 
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